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ABSTRACT 


The principle of virtual work was used to determine 
several typical difference equations governing the dis- 
placements of a taut String subjected To sadist 
transverse load. The integrals for the strain energy and 
the external work were approximated by both the trapezoidal 
mole and Sampson "ss rule. Furvner, the integrand im une 
Perain energy expression was replaced by two different 
mma ce-Oitterence approximations. For ™eemparieon. cwo 
meeneal Githerence equations were cdevelloped irom tie 
Seerrerential equation for the string using central finite 
Peeper oo oes and tne Hermivien ditfterencimae.cneme,. Linc 
finite difference equations, all of which are tridiagonal, 


were evaluated by comparing truncation error. 
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L, JINTRODGerION 


The principle of virtual work is a very usciuleiees 
for finding the displacements of an elastic body produced 
by given forces [ 1] . Im thas thesis Ghevora nears eae 
used to determine several finite difference equations 
meverning the displacements Of a Vault eweiesrre stmt aoe 
moactcrcad FO 2 distributed transverse led¢@ecer Uni tate wear 


mee the strain energy ofesthe string, Usems given» 


H a, \— 
U =3 ; Are dx (i) 


meme PF is the constant internal r1orce inv tie Ssurine. 


fets the displacement, L is the length, and x iS Ghe axial 


@oordinate. The positive directions of F5 ¥> Gy amt & 


7 


mee andicated in Figure 1. 





Pigure 1. A STRING SUBJECTED TO A TRANVERSE LOAD 


AND A TENSILE FORCE 








The pranciple of virtuale cas 
fu = $w Cla2) 
where §U is the variation of U with respect to w and §W 
as the work performed by gq moving through the virtual 
displacement §w. Hence 
&W = [a Swax (Gs) 

ihe trapezoidal rule and Jeanps ens eiuile ss ee emer 
mumericaliy evaluate the continuous antesrais in Houationms 
Hl and 1.3 when the string is divided into kK intervals of 
HMeneth h. The displacements of the strams aresdefined as 
fae Crnd points or nodes of these intervals, and the 
Memiciiueus first derivative in U is replaced by a stakes 
fbpeuierence expression. several approximations to (he wii 
feiawiayive are considered. dhis procedure preduccesa 
Jengtny equation in terms of tne displacements and i1vads 
pomeume Linite number of nodes. Takings the vardacien ole. 
Papo with respect to the displacement at each ete 
nodes’, in accordance wilh the pranciple Of vi reuse om 
Meads to a set of simultaneous linear algebraic equatlenar 
There is one such set of equations for each numerical scheme. 
femepein each set there are one or two typical difference 
equations. The schemes are evaluated by examining the 


mamecation Errors associated with Gach typical equation. 


Itne variation must Saclcivebne CSOMeL lI Ge poumdarny 
Senditions, i.¢€., §w = CG ay x =o and x = L. 








A second procedure for finding Che Grsplecemnen tomes 
an elastic body is to use &@ finite dititerenes appre ine we 
of the governing differential equation, Bien the sa ia 


ae Cd t Pe hen dele CG Wer emma 


aw 
F me + qd = 0 (1.4) 


SevocLtouUbing a finite difference approximacvion for era 
mecond Gaecrivative leads to a set™or samulvancous tiacc. 
emreecpraic equations of tridiagonal form. Two typacal 
mumerec Cliference equavmons are derived=eand Compared. a0™ 
ime bypical equations developed from the virtual work 


Pieineilple. 
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IIT. DEVELOPMENT OF THE FINITE DIFFERENCE EQUATIONS 


USING THE PRINCIPLE OF VIRTUAL WORK 


A. THE TRAPEZOIDAL RULE AND SIMPSON'S RULE 

There are many formulas avaemilablc for numer ico: 
evaluating integrals. Ghe trapezoidal rule and Simpson's 
mle, two of the most commonly usedeeweremchosen sion cai. 
thesis. Both of these find approximations to integrals 
me passing —an interpolating functiom through several modes 
pleeamccorating the interpolating function, ~The trape— 
mordeas Yule uses a straight line as Che interpolating 
myc tion between ie and X44 Where xX, iS eunelva licmen 


S 


Peat the , th MOGe. sonata. ot a7erd as | 2] 


x gees e ; 
Ca b bh @ £00) 
TOK cies 2] 8x59) - £(x5)] — aie on 


Xo 


in which x, < [<x Tne Jas Weim a Se Uchaier eral 


s+1° 
mmciown as the truncation error. Mie inGesral iron. 


wo. x=l, is 
(2o25) 


L 
J meejdx = h Bf (O)HE Ox hee AE CG. 9 HEE) | 
O 
sie Bees 
ae a 
in whichoZ @S4 
im Simpson's rule @ parabola as used as the amverpolau— 


; oe d : 
ing function passing through the nodes ee xO p eile 3s 


igs 





The Aancgegral from x—Oeveesc=—hes 


(2p) 
J, f(x)ax = hb | £(0)448 (oy 428 (4 428 9) 
hoy 
se (x, )+2(2) | ee 


where OS €@ £ L and Kk is even. 


B. DIFFERENCE APPROXIMATIONS TO THE FIRST DERIVATIVE 
The equation containing the conventional central 
muni le Gdalference approximation COn Ene tie weer owen 


at X Seg wae LrUunCcaucvom error is | 2 | 


2/3 
dw = lew ee -h {da-y leg i 
cal a “sti ~ “sa? ~ 2 calm Sete. 


mia whach x <= (sae : 
pe s+] 


Consider the scheme in which Ghe Qilspa@en. wilson. 
marmged act the nodes and Che first derivatives aroma ei 
at the center of the intervals as shown in Figure 2. This 
scheme is referred to as the staggered differencing scheme 
paeeoli-Svuavulon scheme. The central finite et fermen. 
eseuavion for the staggered scheme is derived in Appendix 


fe and is 


* (a>u) 


oe) - 1 (w ~ wi) - h (d=w (2.5) 
— = \"s+1 S 
bi sth h iy) axe Goat 

< < wT 
where Xe > 0544 = Xs41° Hereafter, unless noted otherwise 
Ge denotes that T is located in the interval X_1 bo Seat 


and {544 is in the interval x_ to Xa? and@ is in the 


miceryval x=O CO X=L. 


li 
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(eat cs —— 





We: 


Ke og XS.v2 Ks 54 72 XSer 


Pagure 2. THE STAGGERED” SCHEME 


mnotuner finite difference equation Used Go replaces ea 


merst derivative atv Xo is this higher order expression 


aw) 2 2 ; +h - #70" /arw 
ou) eh eer ee We (are en) 360 2s), 


(2.6) 


Hee corresponding higher order staggered scheme equation is 


+ Laon 
Ch ae -w,) +_h (q ~ Thy fen 
(Be) we Mon 7s) tty Go os) Ee (SO 
STS S+5 
(25) 
Equations 2.6 and 2.7 are derived in Appendix A. 
C, DEVELOPMENT OF THE TYPICAL FINITE DIFFERENCE EQUATIONS 
Mie first typical difference equation is derived usime 
Pema neZOLdal rule aad the cemeral, fanite Gaimerence e-— 
pression. Applying Equations 2.2 and 2.4 to Equation 1.1 


leads to: 


rs 








me 3 2 2 5 (2 
nh a~w eget ax 
#(&) PP roan a2 (te) 
6 ax a Py ae 
(2.8) 


Only those terms in which W, appears are listed in Equation 
2.8 because only the typical equation is of interest. 


faking the variation of Equation 2.8 with respect to the 





displacement We leads to (2.9) 
3 S 
6U . . F (we_p-2w.4w_..)-heF | Ary a" 
Ba re te ees, soils as 
fw P|) onde ax 
S th li Os. sus 


bengne Lhe varlation of the external work with reéespeccrco 


W. leads to 


dW = hq 


§w. : (22 10% 
MNece in Equation 2.9 that the displacements are requared 


at and . This equation does not take into 
s-2 S+2 


consideration the fact that the displacements at Bal and 


a amesevarlable,  Henee. it appears Veaconso le sue. — 
S 


these displacements instead of those at nd 


a e 
S-e s+2 
Consequently, let the interval h be reduced by a factor of 
one half and change the subscripts accordingly to give a 


set of equations that do not omit any nodes. fhe mod11 Le? 


EaiacLons are: 


144 








2 
ST - © (w -2w +w .)-4F a3) 
gW. on Sigal Sk k 1& an 7 ) 
ax3 Zell 
A (2.11) 
fw h 
7 3 ome (2529) 
S 


It can be shown that equations 2.11 and 2.12 are identical 
to the equations developed using the trapezoidal rule and 
the staggered difference scheme. Using the staggered 
scheme eliminates the necessity to halve the internal. 

mS an example of the derivation of the higher Gracr 
PuuMaulons, a typical equation is derived using Simpsons 
molec at eenijiuactien with the higher order staggered scheme. 


Mie ecquation for U based on Simpson's rule can be si ven yin 


the form 
a 2 
Ti eee eee au) +4 (2) oe 
6 obey 6b are 
m (| 
lk Cee 
90 ax (2, 13a) 
Pmiacton 2.13a could also have been siven as 
2 C 
Te re ee +2 (2) +. 
6 ax , cbr ; 
= S+> 
m fh (a) 
~ "p02 SKE 
es 
21 S18) ) 


ls. 








The typical equations associated with Equations 2.13a and 
©.13b are different. The equation derivedwervcm Equawuon 
2.13a is referred to as the 2-4-2 scheme and the equation 
derived from Equation 2.13b is referred to as the 4-2-4 
Ee mMeiie. 


substituting Equation 2.7 into Equation 2.13a gives 


Bust al ; 
ae Peace = (w W 


h 
S s-1) oh (d_ - oa 





Taking the variation of U with respect to w leads to 
S 


i -->. 7 _ A 7 
(W 7 3H ++ eweay) 36 ‘cone aL 2447? 


Ta! OO. Te 2s [f, 5 - 2 (2) 
Y B60 ax /T dx? 1 ee (25159) 


2 





mee Variation of the external work with respec co lbs t Ox 
the integral evaluated by Simpson 's rule with the 2-4-2 


eememe yields 


i: me ove (21s 


le 








All of the typical difference equations that were derived 
using the principle of virtual work are piven in Teblecum 
and II. Table I contains the equations derived from the 
trapezoidal rule and Table II contains the equations derived 
moon olNpsen Ss rule. 

The truncation errors listed in Tables I and II are only 
mart of the total truncation error. Inhese trincaptowmer en. 
eee associated with the finite diitterence evaluabiense en 
derivatives in the strain energy integral. Note that they 
memcisc Of tne ditierence of two derivativesm Yiewcon on 
Peavenee, they cam be converted to a Nigner Orden sots ae 


mexomole, from the first equation of Tablegi 


6 3 
d-w d~w 
tao 7 (a=) co | 
oe ‘€ : ( 
a 
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IITl., THE DIFFERENTIAL EQUATION APPROACH 


A, CENTRAL FINITE DIFFERENCE FORMULATION 

The differential equation for Che Sturpimicetse sm. ae) 
Equation 1.4. Using the conventional central finite difference 
expression to evaluate the second derivative at s leadsero 
[ 2] 


F( ow +0 +) sno + mn ay 
W ra a L 
@aj eu qd. a5 ax 7 (3.1) 


B. HERMITIAN FINITE DIFFERENCE FORMULATION 
Hoe Hermitian formula for the second cderivati vow. 


feniyed in Appendix C, is 


oa = lay é; d 7% ' Gd Ww d —W ' ' 
wow HW =2/( + 10 ( + | 
\ s-l s a? | axe ae) cre | 
Ss-i S Son 
me 6 
_ = 2 ae Cae) 
eho  |ax0 : 
Applying Equation 1.4 at ae eects tO Hquablemeen2 
siail $+] 
gives 
E(w - ew +w cee +10q + 
( S- Ss a 12 (he dk : ae fyag! 
4 6 
ee ( :) = 0 (S83) 
2Ho Vax°/T 


Mameacion 3.3 is the typical dattierence equation for the 


Hermitian finite difference scheme. 


ae 








Equations 3.1 and 3.3, the two typical difference 


éduations; are listed in apie w er. 


out 
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SdUSISTIIG eatutly [e7ruse9 
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IV. COME SHONs 


All of the equations derived give tridiagonal matrices 
centered about the principal diagonal. For a given number 
of nodes, a computer program using these equations should 
have approximately equal storage requirements, computation 
time and round-off error. 

There are two ways to decrease the truncation error. 

One is to use higher order difference equations and the 

Semer 18 to reduce the interval, making more nodes. Nornawa- 
begun of these cause greater round-off error, storage 

me juUEreMments, and computation time, Since the farst amereaces 
the number of non-zero terms in the matrix and the second 
macreases Che size or the matrix. 

The higher order equations derived for this thesis do 
Mmoep 1ncrease the size of the displacement matrix but do 
change the load matrix from a diagonal to a tridiagonal 
Meer xX . 

The object of any numerical analysis method is to obtain 
as accurate an answer as possible within the shortest com- 
Pieacion time. Of the two finite difference solutions to 
the differential equation, the Hermitian method is the better 
erune two from the standpoint of truncation error. 

For the equations derived from the energy solution, 
the truncation error is composed of three parts. Two parts 


are associated with the evaluation of the strain energy and 


23 








work integrals and are of O(h=) for the trapezoidal rule 
and O(n’) for Simpson's rule. The third is the one 
associated with the evaluation of the derivatives in the 
strain energy integral and is Jisted in Tables 1 andy 
This part was of onan when the Central finite dite ewene. 
expression was used and of O(n’) when the higher order 
finite difference expression was used. 

Brom the comparison of the ecquataen. made apis 
meesis by partial truncation error 10 would appeareriaeeree 
pienmoy eduations derived using Simpson's rule andro 
higher-order finite difference equations are the better 
eauacions. However, it must be stressed Chat this evalegaunen 
iemoased Only on Che apparent truncation eCrror anc is. 
Meeomeotmore, NOG Conclusive, 

There is some difficulty in establishing tne total 
Gmeuncation error associated with the equations derived irom 
wee enersy approach. Consequently, it is recommended that 
these equations be evaluated on the computer so that a com- 
plete comparison can be made not only of truncation error 


but also of storage requirements and computation time. 


ou 








APPENDIX A 


THE HIGHER ORDER FINITE DIFFERENCE EQUATION POR THE FIRST 
DERIVATIVE 


Wrens tae Taylor er ito ome eemsy and taking the 


Sac 
mMerivyavive LWiCe fives 


2 
i) = | as | rl »(s) ie i ay) * ef aw | fe 
ax s+] Ax” 7 . ax? S _ ae ee sie axl. 


(A.1) 








end for We gives 


2 2 4 
(3) (5) 1S) AS) 18): 
1 ax | Tebis | ax 2 1dax 6 ax? 


rt 
_ 





(A.2) 


Subtracting Equation A.e from Equation A.1 gives 


fe _[aew - nf ow (es) oes 
ax@ s+] ax@ ae \ ax3 = me abe - 


Gs) 
Solving this —— gives 
see 
Ss 
& 
a>w) 1 (faew amu ne (eae 
oe a, -| 5 - = 4,5 clas aeltetee 
dx ee Vax dx dx 
S S+1 s-l1 S 
(A.4) 


— 








Substituting Equation A.4 into the central finite 
difference equation, Equation 2.4, for the non-staggered 


ater aez tLe lal iets: 

aw bal (w a cae Pan Je 

oo, Stl 12 \ ax® “360: ae an 
S s+1 3 dx 


mise differential equation for the string is 


(ex ee er 
a a vs 
Cs) F (A.6) 


Substituting Equation A.6 into Equation A.5 and truncatinge 


SiVeS 


1 
“ig 
é 


dx on | Stl s-l’ yop "stl “s-1 “Gee (is a 
(A.7) 


moweation A.f is the higher order finite difference equation 


mor the Non-stasesered scheme, 
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APPENDIX B 


THE DERIVATION OF THE FORMULAS AND BQUATIONS USED WITH 
THE STAGGERED SCHEME 


A, THE CENTRAL FINITE DIFFERENCE EQUATION FOR THE FIRST 
DERIVATIVE 


The Taylor series for w Ls 
S+1 


z n! ay) Pak a? —_ 
304 faxt} =, 3840 lax? J, 
STS Sta 


(feeds) 


and Taylor series for w is 
Ss 


2 2 5 3 
7 h as | nh (2 r —h (2 “) 
W ae W —=: - -_ + 
g s+i 2 (2: s+$ 8 Giox ) — 48 dx 3 1 


ne [2s n> ee | 
ate errs” eis zi a ae 5 a5 ® ¢ s 
384 \ax , 3840 \ dx , 
S+ 


2 STZ (Be ca 


poupuracting Mquation B.2 from Equation B.1 gives 


$ 5 ? 
W - W = n=) ae n> & “) {eG (=) T 9) ee 
Stl os Xe ea sth 19eO \ dx st+3 


(B.3) 


eae 








dw 
DO] Valen das 1 Or e gives 


il 
sS+rs 


() = Als (w w.) o ah, (3) £ so) 
d h s+” - > 
Xda ax yd 1920 


(B.4) 
By truncating after the term in w the central finite 
mire rence CQUaLLON Ter Che Dizau deri yvatavc of ome. 


staggered interval is 


where x Se ae ee (Enos) 


eee Bee HTCHER ORDER BPINTIE DLPPERENCE EQUATION FOR FPIRS® 
DERIVATIVE 
[robe ioe one Lavo Series fOr and taking the 


Se 
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Ewoebicvutving Equation B.9 into Equation B.4 sives 


2 
aw) 0 la )-h| fav) _ fee 
dx a LSS Aly | a= axe 
S+3 S+1 = 


60 ED as 
Dil Xx SS 
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Equation B.12 is the higher order finite difference 
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equation for the staggered scheme. 
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APPENDIX C 
THE DERIVATION OF THE HERMITIAN FINITE DIFPERENCK EQUATION 
FOR THR SECOND DERIVATIVE 


Adding Equations A.1 and A.2@ gives 
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The Taylor series derived expression for the second 
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Substituting Equation C.2 into Equation C.3 gives 
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